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Abstract. In the Poynting Flux-dominated outflow (the initial ratio of the electromagnetic energy flux to the 
particle energy flux uq 3> 1) model for Gamma-ray bursts, particularly the 7— ray emission phase, nearly half of 
the internally dissipated magnetic energy is converted into the 7— ray energy emission and the rest is converted 
into the kinetic energy of the outflow. Consequently, at the end of the 7— ray burst, a decreases significantly 
(a ~ 1 or even smaller). We numerically investigate the very early reverse shock emission powered by such mildly 
magnetized outflows interacting with medium — uniform interstellar medium (ISM) or stellar wind (WIND). We 
show that for a ~ 0.05 — 1 and typical parameters of Gamma-ray bursts, both the ISM-ejecta interaction and the 
' WIND-ejecta interaction can power very strong optical emission (i«h ~ 10 — 12th magnitude or even brighter). 

Q\ . Similar to the very early afterglow powered by the non-magnetized ejecta interacting with the external medium, 

C*") 1 the main difference between the ISM-ejecta interaction case and the WIND-ejecta interaction case is that, before 

the reverse shock crosses the ejecta, the R-band emission flux increases rapidly for the former, but for the latter 
, it increases only slightly. 

At the very early stage, the ejecta are ultra-relativistic. Due to the beaming effect, the random magnetic field 
generated in shocks contained in the viewing area is axisymmetric, unless the line of sight is very near the edge 
of ejecta. The formula il nct ~ 0.606 2 /(l + b 2 ) (where b is the ratio of the ordered magnetic field strength to that 
|' of random one) has been proposed to describe the net linear polarization of the synchrotron radiation coming 

1~ \ from the viewing area. For a ~ 0.05 — 1, the ordered magnetic field dominates over the random one generated 

in the reverse shock (As usual, we assume that a fraction eb ~ 0.01 of the thermal energy of the reverse shock 
has been converted into the magnetic energy), the high linear polarization is expected. We suggest that the linear 
polarization detection of the early multi-wavelength afterglow is required to see whether the outflows powering 
GRBs are magnetized or not. 
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1. Introduction been various implementations of this concept (Usov 1992, 

1994; Thompson 1994; Blackman, Yi & Field 1996; Katz 
It is widely accepted that 7 -ray bursts (GRBs) are pow- 19g7; Meszaros & Rees 1997) . The Poynting flux model 
ered by the dissipation of energy m a highly relativistic ig ^ Ught of ^ foUowing facts (see zhang & Meszaros 
wind, driven by gravitational collapse of a massive star (20Q4) for & recent review) . The Poynting flux outflow 
into a neutron star or a black hole (see Meszaros 2002 for can transport a large amount of energy without carrying 
a recent review). As the observed emission is powered at a many baryons It can a n ev iate the inefficiency problem of 
distance far from the central source, key questions remain the intemal shock model and can alsQ alleviate the mag _ 
unanswered. One of them is the Gamma-ray burst engines netic fieJd amplification prob i em in GRBs and afterglows, 
(see Cheng & Lu 2001 for a review). In the standard fire- n provides the possibility of achieving narrow "peak en- 
ball model, the Gamma-ray burst are powered by the col- ergy „ distributions . In part i C ular, the Poynting flux model 
hsions of baryon dominated shells with variable Lorentz provides the most natural explanation so far for the very 
factors (Paczyhski & Xu 1994; Rees & Meszaros 1994). high polarization during the 7 _ ray emission phase of 
However, Poynting flux-driven outflows from magnetized Qpjg 02 i 2 06 (e.g. Coburn & Boggs 2003; Lyutikov, Pariev 
rotators is another plausible explanation and there have & Blandford 2003; Granot 20 03; However, see Rutledge 

~ ; ~ : „ _ & Fox (2004) for argument), although other alternative 
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explanations remain (e.g. Shaviv & Dar 1995; Waxman 
2003). 

In the past several years, many publications have fo- 
cused on the dissipation (the energetic non-thermal 7— ray 
emission) as well as the acceleration of the Poynting 
flux outflow (e.g. Usov 1994; Thompson 1994; Smolsky & 
Usov 1996; Meszaros & Rees 1997; Lyutikov & Blackman 
2001; Spruit, Daigne & Drenkhahn 2001; Drenkhahn 2002; 
Drcnkhahn & Spruit 2002). In this paper we turn to in- 
vestigate the very early reverse shock emission powered by 
such a magnetized outflow, just as Sari & Piran (1999) and 
Meszaros & Rees (1999) have done for the baryon domi- 
nated fireball. One thing inciting us to do this is that mod- 
cling the very early afterglow of GRB 990123 and GRB 
021211 suggests the reverse shock emission region is mag- 
netized (Fan et al. 2002; Zhang, Kobayashi & Meszaros 
2003). 

In the "internal" magnetic dissipation model, GRBs 
are powered by the magnetic energy dissipation at a ra- 
dius r ~ 10 13 — 10 14 cm. At the end of the 7— ray burst, a 
significant fraction of magnetic energy has been dissipated 
by magnetic reconnection or other processes. Nearly half 
of the dissipated magnetic energy has been converted into 
the 7— ray emission and the rest has been converted into 
the kinetic energy of the outflow (see Spruit & Drenkhahn 
2003 for a recent review). Consequently, a decreases sig- 
nificantly (<~ 1 or even smaller). At a much larger radius, 
where the outflow begins to be decelerated significantly, 
the reverse shock emission (the very early afterglow) is 
expected; this is what we focus on. 

At the final stage of the preparation of this manuscript, 
a paper by Zhang & Kobayashi (2004) appeared. In that 
paper, the very early reverse shock emission from an arbi- 
trary magnetized ejecta has been analytically investigated. 



2. The mild-magnetized outflow 

As mentioned before, there are lots of publications focused 
on the acceleration of the magnetized outflow. One of them 
is Drenkhahn (2002), in which part of the magnetic en- 
ergy coupled with the outflow is dissipated internally by 
reconnection and the Lorentz factor of the flow increases 
steadily with radius (r oc r 1 / 3 ). Here we do not discuss 
that topic further and just take the numerical example 
presented in Drenkhahn & Spruit (2002) as the starting 
point of our calculation: at the end of the prompt 7— ray 
emission phase, the bulk Lorentz factor of the outflow is 
7] ~ 300; the ratio of the electromagnetic energy flux to 
the particle energy flux, a ~ 0.05 — 1 (In principle, much 
lower a is possible, for which the reverse shock emission 
is similar to that of the usual fireball, which is beyond 
our interest); the total kinetic energy (including the mag- 
netic energy) is of the order of the typical 7— ray emission 
energy, i.e., E kin ~ 10 53 ergs. 



3. The reverse shock emission 

3.1. The dynamical evolution of ejecta 

Generally, the dynamical evolution of the ejecta can be di- 
vided into two phases — (i) Before the reverse shock crosses 
the ejecta, i.e., R < R cro (R is the radial coordinate in 
the burster frame; i? cro is the radius at which the reverse 
shock crosses the ejecta); at that time two shocks exist. 
The dynamical evolution of the ejecta is governed by the 
jump condition of shocks (e.g., Blandford & McKee 1976; 
Sari & Piran 1995); (ii) After the reverse shock crosses the 
ejecta, i.e., R > i? C ro, in this case only the forward shock 
exists. The hydrodynamical evolution can be calculated 
by taking the generic dynamical model of GRB remnants 
(e.g., Huang, Dai & Lu 1999; Huang ct al. 2000; Feng et 
al. 2002). 



3.1.1. The dynamical evolution for R < R CTO 

Similar to Sari & Piran (1995), the dynamical evolution 
of the ejecta is obtained by solving the jump condition for 
strong shocks. For the ejecta interacting with the external 
medium, there are two shocks formed, one is the forward 
shock expanding into the medium, the other is the reverse 
shock penetrating the ejecta. There are four regions in 
this system: (1) The un-shocked medium; (2) The shocked 
medium; (3) The shocked ejecta material; (4) The un- 
shocked ejecta material. The medium is at rest relative to 
the observer. The bulk Lorentz factors Tj (j — 1,4, T4 = 
77) and the corresponding velocities 0r } = (1 — 1/T?) 1 / 2 
are measured by the observer. Thermodynamic quantities: 
rij, pj, ej, B- (particle number density, pressure, internal 
energy density, magnetic field strength) are measured in 
the fluids' rest frame (we assume the un-shocked ejecta 
and medium are cold, i.e., e\ = e\ = 0), so is the pe,i (the 
magnetic pressure). The equations governing the forward 
shock are (Blandford & Mackee 1976) 



n2/n 1 — 4r 2 + 3, e 2 /n 2 = (T 2 — l)m p c 2 



(1) 



Below, following Kennel & Coroniti (1984, hereafter 
KC84) we try to derive the (90°) shock jump condition 
governing the reverse shock with the MHD conservation 
laws by assuming the magnetic field frozen in the outflow 
is nearly toroidal (KC84) : 



n 4 u 4 = n 3 u 3 , 



74M4 + 



EB4 



Pa 

H4U4 H h 



73M3 + 



TI4U4 8TTT14U4 



EB3 

4irri4U4 ' 

^3 

H h 



U4U4 87rn4U4 ' 



(2) 
(3) 



(4) 



where _B; and E denote the shock frame magnetic and 
electric fields respectively, E = 7J4B4/74 = u 3 B 3 /j 3 . 71 
(i = 3, 4) is the Lorentz factor of the fluid measured in 
the reverse shock frame, uf = 7? — 1, $ = 7*1/7;. I 1 is 
the specific enthalpy, which is defined by m = m v c 2 + 
7iPi/[(7i — l)rii] for a gas with an adiabatic index 7;. 



Y. Z. Fan, D. M. Wei and C. F. Wang: Reverse shock emission powered by magnetized ejecta 



3 



Solving equation (3) for p 3 and inserting the resulting 
expression into equation (4) leads to 



a(Y 2 - 1) + 



2[1 + ct(1-F)] Ui 



73 - In 



M373 74 73 

2P 4 u 4 2(73 - l)m p c 2 



74 



= 0, 



(5) 



n 4 ^47f u 3 73M47I 

where a = B\j '[An^^-j 2 ] and Y = With equation 

(3), the downstream pressure P3 can be calculated as fol- 
lows: 



P3 = 



73 - l r 74 



73 73 



\n-i[L4- 



(6) 



For 74 3> 1 and P4 — 0, equation (5) can be arranged 
into (Assuming 2(73 — IVU37374 can be ignored) 



(l + q)(2- 73 ) 4 , , 73 -4 2 t| -273 + 2 

"3 + I^T C To cr 



73 
,73 - 1 



4 73 

)>3 + (%-^) 2 - 2 



7l 



0. 



(7) 



73 ' ' " 273 
For 73 =4/3, the above equation can be greatly simplified 

8(1 + o)u\ - (8cr 2 + 10cr + l)u\ + cr 2 = 0, (8) 

whose solution is equation (4.11) of KC84. 

The total pressure in region 3 can be calculated by 



P 



3. tot 



p 3 



2 R /2 

4-° " 



. The equality of pressure and veloci- 



ties along the contact discontinuity yields P3,tot = P2 = 
4Flnim p c 2 /3, T 2 = T 3 . For the Lorentz factor of the 
reverse shock (measured by the observer) r rs h 1, 74 
can be expressed as 74 w (?7/r rs h + T TS i 1 /r])/2, which in 
turn yields r rs h ~ (74 — 114)7]. On the other hand, 73 
can be expressed as 73 ss (r2/r rs h + r rs h/r2)/2, which 
in turn yields T2 ~ (73 + M3)r rs h, where the relation 
T 2 = T3 has been taken. Combing these relations we have 
1^2 ~ (73 + U3X74 — U4)v- Finally, we have the equation 
(the equality of pressure) 



B'l u 4 (73 - 1) 74 



47T(7 ^ U 3 73 



[^(l + a(l-Y))- 

73 



TOpC 



J + 2u 2 1 



M4 ^"3 

4[(73 + "3)(74 - U4)rf\ 2 nim v c 2 



(9) 



where B 4 = 2.6 x 10 20 G R- 1 ^ 1 . Equations (1), (5) and 
(9) are our basic formulae, with which we can calculate the 
dynamical evolution of the magnetized outflow interacting 
with the medium, and then the reverse shock emission. 
However, these equations cannot be solved analytically 
unless a ^> 1. For a ~ 1 or smaller, only the numerical 
calculation can be performed, which is to be presented at 
the end of this section. 

R cro can be determined as follows: /3 rs h) the velocity of 
the reverse shock in the observer's frame, can be parame- 
terized as (Sari & Piran 1995) 



A-sh — 



r 3 n 3 - r 4 n 4 



(10) 



Differentially, (/3r 4 — (3 rs h)dR = dA, where A is the 
width of the reverse shock penetrating into the ejecta 
measured by the observer. R cro is determined by J dA = 

J Rcro (/3r 4 — Prsh)dR — A , where A is the width of the 
ejecta measured by the observer, which can be estimated 
as A w cTg /(l + z) (Tg is the observed duration of 
GRBs). 

3.1.2. The dynamical evolution for R > R cro 

After the reverse shock has crossed the ejecta, only the 
forward shock exists, whose dynamics have been discussed 
in great detail (e.g. Huang et al. 1999; Huang et al. 2000; 
Panaitescu & Kumar 2001; Feng et al. 2002). However, 
in the current work, the ejecta is magnetized and how to 
convent the magnetic energy into kinetic energy is poorly 
known. But the energy conservation must be satisfied. As 
a zeroth order approximation, here we take Huang et al.'s 
(1999) differential equation to depict the dynamical evo- 
lution of the magnetized ejecta 



dT_ 

dm 



M'- + em + 2(l-c)Tm' 



(11) 



where T is the bulk Lorentz factor of the outflow, (3 — 
(1 — 1/r 2 ) 1 / 2 is the corresponding velocity in unit of light 
speed; m is the mass of the medium swept by the ejecta, 
dm = 4:irR 2 nim p dR; < e < 1 is the radiation effi- 
ciency; M' c - } is determined by (the energy conservation) 
7croA^ojC 2 = (1 — e ra d)Poj where e ra( j is the fraction of to- 
tal energy radiated for R < i? cro , Eq is the initial isotropic 
energy of the ejecta. 7 cro is the bulk Lorentz factor of the 
ejecta at i? C ro, which can be determined by equations (1), 
(5) and (9). 

3.2. The synchrotron radiation 
3.2.1. Electron distribution 

In the absence of radiation loss, the distribution of the 
shock accelerated electrons behind the blastwave is usually 
assumed to be a power law function of electron energy, i.e., 



dK 
djc 



OC 7 C P , (7e,m < 7c < 7c,m), 



(12) 



where p w 2.2, 7 CjM = 10 8 (B'/1 G ) _1/2 is the maximum 
Lorentz factor (Dai, Huang & Lu 1999), where B' is the 
comoving frame magnetic field strength. As usual, for ma- 
terial heated by the forward shock, we assume the mag- 
netic energy density in the comoving frame is a fraction 
£b of the total thermal energy, B' 2 can be estimated as 



B 



= e B (47 2 + 3)(7 2 - l)mm p c 2 



(13) 



where 7 2 = T 2 for R < R cro and 7 2 = T for R > R cro . For 
region 3 (R < R ao ), B' 3 = u 4 B 4 /u 3 . 

Similarly, for region 2 and 3, the thermal energy of 
electrons in the comoving frame is assumed to be a fraction 
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e e of the total thermal energy, then 7 C m can be estimated 
by 1 



P - 2 m p , 



(Y-l). 



(14) 



p - 1 m e 

where m c is the rest mass of electron, 7' = 72 for region 2 
and 7' = 7 P: th + 1 for region 3, where the thermal Lorentz 
factor of the shocked proton 7 p , th can be estimated by (see 
equation (6)) 



7 P ,th = e 3 /n 3 m p c 2 



1 ^-(i + .(i-y))-^]^Mi5) 



73 73 M4 m p C 

It is well known that radiation loss may play an impor- 
tant role in the process. Sari, Piran & Narayan (1998) have 
derived an equation for the critical electron Lorentz fac- 
tor, 7c, k (k = 2, 3 represent region 2 and 3), above which 
synchrotron radiation is significant 

127rm c r'(l - (3'n)c 



7c,k 



(16) 



where t is the observer time, fi = cos 6, is the angle 
between the velocity of emitting material and the line of 
sight; Throughout the rest of this work, T' = T 2 (/?' is the 
corresponding velocity in unit of light) in the presence of 
reverse shock and r" = L at later time. If the emitting 
material moves on the line of sight, equation (16) reduces 
to the familiar form 7 C k = 6 VPSri t ■ For electrons with 

' ' (TtI r>^ t 

Lorentz factors below 7 c .k, the synchrotron radiation is 
ineffective. For electrons above 7 c ,k, they are highly radia- 
tive. 

In the presence of steady injection of electrons accel- 
erated by the shock, the distribution of electrons with 
7c > 7c, k has a power law function with an index of 
p+1 (Rybicki & Lightman 1979), while the distribution of 
adiabatic electrons is unchanged. The actual distribution 
should be given according to the following cases: 

(i) For 7 C < 7c, m < 7 c ,m, i-e., the fast cooling phase 

(7c < 7c < 7e,m), 



dK _ r J 7c~ 2 , 
dr7 " 1 1 



{ " >r ': (7o,m < 7c < 7c,m), ^ 



r ,1 P-I 1 7£ V ,.!,, 
^1 = 1 J JV tot j 



7c 



P 7c 



P 



(18) 



where N tot is the total number of radiating electrons 
involved. 

(ii) For 7c, m < 7 C < 7c.m> i- e -> the slow cooling phase 



d> I 70 (P+1) 



(7e,m < 7c < 7c), 

(7c < 7c < 7c, m), 



where 



C 2 = [ 



7, 



i-p 



7 i-p 

/e,m 



p-i p(i-p) 



7c7 e ,M 
P 



(19) 



(20) 



1 Here the possible but poorly constrained e pair generation 
during the 7— ray emission phase and its impact on the reverse 
shock emission (e.g. Pilla & Loeb 1998; Li et al. 2003; Fan & 
Wei 2004; Fan et al. 2004a) have not been taken into account. 



As the reverse shock has crossed the ejecta, there are 
no freshly heated electrons injected and the ejecta cools 
adiabatically. Here we investigate the decay of the mag- 
netic field and the cooling behavior of electrons. In the 
current case, the magnetic field is ordered and satisfies 
the magnetic flux conservation. As usual, the ejecta are 
in the spreading phase and the width can be estimated 
as A' - R/T. Therefore, B' 3 oc R^T^/A' oc R~ 2 . The 
cooling behavior of electrons is more difficult to estimate. 
Here we simply assume the cooling of electrons is not much 
different from that of the non-magnetized case. According 
to Meszaros & Rees (1999), 7 c , m oc R-( 9 - s )/ e ( so does j c , 
since all of them cool by adiabatic expansion only), where 
s = for the ISM and s = 2 for the stellar wind. In the 
case of fast cooling, once ^ c . b s = e B' 3 r ) 2 T 3 /2Trm c c drops 
below the observed frequency, the flux drops exponentially 
with time (Sari & Piran 1999). If the equal arriving time 
surface has been taken into account, the flux drops more 
slowly. 

3.2.2. Relativistic transformations 

In the co-moving frame, synchrotron radiation power 
at frequency v 1 from electrons is given by (Rybicki & 
Lightman 1979) 



7e 



(21) 



where e is electron charge, v' c — 3-y 2 eB' / \Airm e c) , j' c 
min{7 c , m ,7c} and 



r+00 

F{x)=x K 5/3 (k)dk, 

J X 



(22) 



with K 5 / 3 (k) being the Bessel function. We assume that 
this power is radiated isotropically in the comoving frame, 

dP'(u') _ P'{v') 

The angular distribution of power in the observer's 
frame is (Rybicki & Lightman 1979; see also Huang et 
al. 2000) 



1 



dP'{v') 



1 



P'{v') 



dP{v) _ 

dn ~ L' 3 (l-/?» 3 dSV ~ r /3 (l-/3'/x) 3 4tt 



(23) 



(1 + z)T'(l - (3'ii) 



(24) 



where z is the redshift of the ejecta. Then the observed 
flux density at frequency v is 



S v = 



1 (dP{v) A 
A V dfl Dl 



T' 3 (l - (j'nY 4tt£> 



P'((l + z)T'(l-f3'^), (25) 



where A is the area of our detector and Dl is the lu- 
minosity distance (we assume Ho = 65km s _1 Mpc -1 , 
Vt M = 0.3, Q A = 0.7). 



Y. Z. Fan, D. M. Wei and C. F. Wang: Reverse shock emission powered by magnetized ejecta 



5 



3.2.3. Equal arrival time surfaces 

Photons received by the detector at a particular time t are 
not emitted simultaneously in the burster frame. In order 
to calculate observed flux densities, we should integrate 
over the equal arrival time surface determined by (e.g. 
Huang et al. 2000) 



' = / — J? — dR = const, 



fVc ' (26) 

within the boundaries. 

3.3. Numerical results 

In our calculation, the number density of the medium (in 
unit of cm -3 ) has been taken as 



const. , 

3.0 x 10 35 A,i? 



(ISM), 
(WIND), 



5 M yr- 1 V 10 a kn 



(27) 



rr) -1 , M is 



respectively, where A* = -^p 

the mass loss rate of the progenitor, w w is the wind velocity 
(Dai & Lu 1998; Chevalier & Li 2000). In our numerical 
calculation, we take A* = 1. 



For illustration, we take 



10 bd ergs, 



1, 



L = 2 x 10 51 ergs s" 1 (equally, T 90 = 50(1 + z)s), r) = 300, 
p = 2.2, P 4 = and e c = 0.3. In the case of ISM- 
ejecta interaction, the protons in the region 3 are only 
mild-relativistic or even sub-relativistic, but electrons are 
ultra-relativistic, so 73 = 13/9. In the case of WIND-ejecta 
interaction, protons heated by the reverse shock are rela- 
tivistic, so 73 = 4/3. 

3.3.1. ISM-ejecta interaction case 

The sample very early R-band (i/r = 4.6 x 10 14 Hz) light 
curves have been shown in figure 1. Before the reverse 
shock crosses the ejecta, electrons accelerated by the re- 
verse shock are in the slow cooling phase and the syn- 
chrotron emission at R-band increases rapidly with time. 
At tens of seconds after the main burst, the reverse shock 
emission at R band is bright to wr ~ 10 — 12th magni- 
tude and the forward shock emission is relatively dimmer. 
Therefore, the very early reverse shock emission can be de- 
tected independently. After the reverse shock has crossed 
the ejecta, there are no freshly accelerated electrons in- 
jected and the R-band emission drops sharply. 

In figure 1, there are two interesting phenomena: (i) 
The peak flux at R-band increases with the increasing a 
for a < 0.1, but for a > 0.1 the peak flux at R-band de- 
creases with the increasing a (A similar result has been 
obtained by Zhang & Kobayashi 2004). This behavior 
can be understood as follows: For a <C 1, the electrons 
heated by the reverse shock is in slow cooling phase. At 
i? cro , the typical synchrotron radiation frequency v m , hs — 
T27c m e-B 3 /[27r(l + z)m c c] is much lower than j/r, as is 

t'c.obs- We have approximately F VR cx B' 3 ^ 1 '' 2 . Roughly 
speaking, the observed flux increases with the increasing 




t(8) 



Fig. 1. The very early R-band (vr = 4.6 x 10 14 Hz) light curve 
powered by the mildly magnetized outflow (the degrees of the 
magnetization have been marked in the figure) interacting with 
the interstellar medium. The parameters taken in the calcula- 
tion are: 2 = 1, E kin = 10 53 ergs, L = 2 x 10 51 ergs s"\ rj = 300, 
ni = 1cm -3 , e c = 0.3 and the radiation efficiency e = e c . For 
C7 = and the forward shock, it is assumed that ea = 0.01. 



a. For larger a, the reverse shock has been suppressed 
and the electrons involved in the emission decrease. So 
F VR drops again (see Zhang & Kobayashi 2004 for more 
detailed explanation), (ii) For a ~ 1, the crossing time 
t cro (at which the reverse shock crosses the ejecta; in fig- 
ures 1 and 2, it equals the peak time of the reverse shock 
emission) is much shorter than T90. Here we explain this in 
some detail. In the presence of reverse shock, differentially, 
fcro satisfies 2 



dtcYO (I 



1-0T 



dA 



with r rsh » (73 

Approximately 

^cro / dt ci 



~ Prsh C 

it3)r3, we have dt. 



(1 + z)- rsh 



F 2 
1 3 



1+z 



(28) 



1 



(73,cro + U3,cro) 2 



0(1 ■ 



(J3+U3) 2 c 



(29) 



where 73,01-0 and 1*3, C ro are the corresponding value of 73 
and M3 at R cm . In the case of non-magnetization and the 
reverse shock is relativistic, 73.0-0 ~ y/9/8, which yields 
^cros ~ (1 + z)Ao/2c sa Tgo/2 (This result coincides with 
the "classical" result of Sari & Piran 1995). While for 
<T = 1, 73. C ro ~ 1-37 which results in i cros w Tgo/5, which 
is quite consistent with our numerical result in Fig.l. 



2 If the reverse shock is Newtonian, i.e. 
where < / < 1. We have /3r 4 - Prsh « 

should be written into dt CTO 



r rs h 
X 



(1 - /)r 4 , 

Equation (28) 



2/ 



Approximately, t CIO 



4jf » T 90 . 
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Fig. 2. The very early R-band light curve powered by the 
mildly magnetized outflow (the degree of the magnetization has 
been marked in the figure) interacting with the stellar wind. 
The parameters taken here are the same to those of figure 1 
except m = 3 x 10 35 i?~ 2 cm~ 3 . 

Equation (29) applies to the following WIND-ejecta inter- 
action case as well, if only the reverse shock is relativistic 
or at least mild-relativistic. 



X 




Fig. 3. Coordinates used in the calculation of the polarization 
properties of a slab of mixed field (after Laing 1980). 

However, now F UR oc B' 3 1 . Consequently, the reverse 
shock emission powered by the high a ejecta interacting 
with WIND is dimmer than that powered by the low a 
ones. 



3.3.2. WIND-ejecta interaction case 

In the case of WIND-ejecta interaction (see figure 2), be- 
fore the reverse shock crosses the ejecta, the electrons ac- 
celerated by the reverse shock are in fast cooling phase 
and the R-band emission increases only slightly with time. 
This temporal behavior is very similar to that of non- 
magnetized fireball case (See Wu et al. (2003) for an ana- 
lytical investigation). At t cro , the reverse shock emission at 
R band is very bright (tor ~ 9— 10th magnitude) and the 
forward shock emission is relatively dimmer. Therefore, 
the very early reverse shock emission can be detected in- 
dependently, too. After the reverse shock has crossed the 
ejecta, the R-band emission drops sharply. 

Since the current reverse shock is relativistic, as im- 
plied by equation (29), t CTO decreases with increasing a. 
However, in figure 2, the R-band reverse shock emission is 
brightest at a = 0, which seems to be inconsonant with the 
result shown in figure 1. The main reason for this "diver- 
gence" is: For R < 10 17 cm, the WIND is far denser than 
the ISM. Consequently, the reverse shock is very strong 
and the ejecta has been decelerated significantly at a ra- 
dius ~ 3 x 10 15 cm (Note that in the case of ISM-ejecta 
interaction, the corresponding radius is ~ 10 17 cm). Even 
for a = 0, the electrons heated by the reverse shock are 
in the fast cooling phase and ^ mj0 bs is much higher than 
the observer frequency i>r. With the increasing er, the 
corresponding R clo decreases. As a result, B' 3 increases. 



4. The linear polarization 

It is well known that for the totally ordered magnetic con- 
figuration, high linear polarization is expected. For the 
random magnetic configuration, mild linear polarization 
can be expected if some specific geometry effects have 
been taken into account. Sometimes the magnetic field is 
not only ordered or only random. Then it is interesting to 
investigate its linear polarization. Here we propose a sim- 
ple formula to describe the linear polarization properties 
of a slab of such a mixed magnetic field, with which we 
can see the impact of the ordered field on the polarization. 

Following Laing (1980; See his Appendix Al for detail), 
the coordinates involved are defined as follows ( see figure 
3): 

a is the angle between the plane of the ejecta and 
the line of sight; x, y, z are rectangular coordinates with 
the z-axis pointing towards the observer (i.e., the direc- 
tion n) and the y-axis parallel to the "local" plane of the 
ejecta; x', y' are coordinates in the plane of the ejecta, 
y 1 is parallel to y; 6 is the angle between the field di- 
rection and the x' axis at any point in the ejecta; x is 
the position angle of the E-vector of the polarized radia- 
tion, measured from the x — y plane. Therefore the ran- 
dom (ordered) magnetic-field vector B ran (B or d) at a point 
in the slab are 3 B ran = i3(cosf?sina, sin (9, cos 8 cos a), 

3 According to Medvedev & Loeb (1999), the configuration 
of the magnetic field generated in shocks is tangled with the 
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B or d = Bq (cos #o sin a, sin6>o, cos #0 cos a) respectively. 
Thus the total magnetic field is 

B = B[(cos8 + 6cos6> ) sin a, sin# + 6sin6* , 

(cos 8 + 6 cos 6*o ) cos a], (30) 

where b = B /B. The electric field of a linearly polarized 
electromagnetic wave is directed along the vector 

e = n x B 

= [—(sin 9 + b sin 9 ), (cos 9 + b cos 9 ) sin a, 0]. (31) 
Now x satisfies 

sin a(cos 9 + b cos 9 ) 



tanx = -- 



sin 6 + b sin O 



(32) 



With equation (32), it is easy to get expressions for cos(2x) 
and sin(2%). With equations (l)-(3) of Laing (1980) and 
assuming the spectra index p = 3, we have the Stokes 
parameters 



Q = -[cos 2 a + 26 2 (sin 2 8 - sin 2 a cos 2 9 )}ttB 2 



U = 1.5sinasin26» 7r6 2 B 2 



(33) 



7TU D , (34) 

I = (l + sin 2 Q! ) + 2fe 2 (sin 2 6' + sin 2 acos 2 6' )7rB 2 . (35) 

For 9 = tt/2 or 37r/2, we have U = and the degree of 
the linear polarization 



IT : 



Q 3 cos 2 a + 2fe 2 3 b 



> 



I 4(l + sin 2 a) + 26 2 " 41 + & 2 ' 



(36) 



For i > 1, II has the maximum value LT — |. In the 
current work, a <~ 0.05 — 1, the corresponding toroidal 
magnetic field is far stronger than that generated in re- 
verse shock, i.e., b 3> 1, so the local point polarization can 
be as high as 75%. For ultra-relativistic ejecta, due to the 
beaming effect, only the emission coming from a very tight 
cone around the line of sight can be detected. If the line of 
sight is slightly off the symmetric axis of the ordered mag- 
netic field, the orientation of the viewed magnetic field is 
nearly the same. The high linear net polarization is ex- 
pected since the local high linear polarization cannot be 
averaged effectively. The detailed numerical calculation of 
the net polarization will be presented elsewhere. 

Here, for simplicity, following the treatment of Granot 
& Konigl (2003), the net Stokes parameters of the or- 
dered magnetic field (ford, Qord, ^ord) and of the random 
magnetic field (U ran , Q 
Therefore 



ran, ^ran) are calculated separately. 



+ <9ord + AT ord nrn b 2 

n„ct = — — — — = - — — — = 0.60 



^ran H~ -^ord /ran H~ ^ord 



i + b 2: 



(37) 



front of shock surface. As a result, locally, in a finite scale, it 
is planar. 



where <5 ra n = U lan = J7 rd = for the symmetrical view- 
ing area, and I ex B 2 for p = 3. For b — > oo 4 , n not ~ 0.60, 
so we take A = 0.60. 

Equation (37) is favored by the fact that for b = 1, 
1/V3 and 0, it gives II nct = 0.30, 0.15 and respectively, 
which coincides with the result of Granot & Konigl (2003) 
excellently. Then we believe that equation (37) provides 
us a rough but reliable estimation on the impact of the 
ordered magnetic field on the linear polarization. Equation 
(37) is valid only when the viewed emitting region for the 
random magnetic field is axisymmetric. If it is not, the 
random magnetic field may play an important role. The 
detailed calculation for that case is beyond the scope of 
this paper. 

5. Discussion and conclusion 

The reverse shock emission in the framework of the stan- 
dard fireball model of GRBs has been discussed in great 
detail (e.g. Sari & Piran 1999; Meszaros & Rees 1999; 
Wang et al. 2000; Kobayashi 2000; Wu et al. 2003; Zhang, 
Kobayashi & Meszaros 2003; Nakar & Piran 2004). The 
very early afterglow of the X-ray Flashes has been investi- 
gated by Fan et al. (2004b) recently. For typical parame- 
ters and reasonable assumptions about the velocity of the 
source expansion, a strong optical flash niR w 9 — 17th 
magnitude is expected (e.g. Sari & Piran 1999; Wu et 
al. 2003; Fan et al. 2004b). However, despite intensive ef- 
forts, only three candidates (GRB 990123, GRB 021004 
and GRB 021211) have been reported (Sari & Piran 
1999; Kobayashi & Zhang 2003; Wei 2003 and references 
listed therein). It is unclear why. Interestingly, model- 
ing the reverse shock emission of GRB 990123 and GRB 
021211 suggests that the reverse shock emission region is 
magnetized — In other words, the magnetic energy density 
in region 3 is far stronger than that in region 2 (Fan et 
al. 2002; Zhang, Kobayashi & Meszaros 2003). There are 
two possible explanations: One is that the magnetic field 
coming from the central source has been dissipated signif- 
icantly, i.e., the case considered in this paper. The other is 
that the magnetic field is generated in internal shock. In 
the internal shock model, the generated magnetic field can 
be as high as 0.01 — 0.1 times the total thermal energy of 
the shocked baryons. The generated magnetic field is ran- 
domly oriented in space, but always lies in the plane of the 
shock front, for which the jump condition derived in §3 is 
satisfied in the coherence scale. More importantly, the an- 
nihilation timescale for the random magnetic field is much 
longer than the dynamical timescale of the fireball, then 
the existing of the generated magnetic field can affect the 
very early afterglow (Medvedev & Loeb 1999). However, 
the coherence scale of the generated magnetic field is so 

4 In Lyutikov, Pariev & Blandford (2003), the Stokes pa- 
rameters are pulse-integrated and the resulting il nct = 0.56. 
Considered that the photons emitted at the same time but 
with different angles arrive at different time, a more detailed 
calculation suggests n net = 0.64. Such a difference is not large. 
Therefore, and partly for simplicity, we take n nct = 0.60. 
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small (~ 10 3 cm) that there is no net polarization in the 
early multi- wavelength emission unless some geometry ef- 
fects have been taken into account (e.g. Medvedev & Loeb 
1999). However, as shown in §4, if part of the magnetic 
field is ordered, high linear polarization can be detected 
(see also Granot & Konigl 2003). Therefore polarization 
detection at very early times may provide us the chance to 
distinguish between the usual baryon-rich fireball model 
and the Poynting flux-dominated outflow model for GRBs. 

The predicted very early afterglow in the R band is 
bright to tor, ~ 10 — 12th magnitude, which is strong 
enough to be detected by current telescopes, such as the 
ROTSE-IIIa telescope system, which is a 0.45-m robotic 
reflecting telescope and managed by a fully-automated 
system of interacting daemons within a Linux environ- 
ment. The telescope has an f-ratio of 1.9, yielding a field of 
view of 1.8x1.8 degrees. The control system is connected 
via a TCP /IP socket to the Gamma-ray Burst Coordinate 
Network (GCN), which can respond to GRB alerts fast 
enough (< 10s). ROSTE-IIIa can reach 17th magnitude 
in a 5-s exposure, 17.5 in 20-s exposure (see Smith et 
al. 2003 for details). Another important instrument for 
detecting the very early afterglow is the Ultraviolet and 
Optical Telescope (UVOT) on board the Swift Satellite. 
The Burst Alert Telescope (BAT) is another important 
telescope, with which hundreds of bursts per year to bet- 
ter than 4 arc minutes location accuracy will be observed. 
Using this prompt burst location information, Swift can 
slew quickly to point the on-board UVOT at the burst 
for continued afterglow studies. The spacecraft's 20—70 
second time-to-target means that about ~100 GRBs per 
year (about 1 /3 of the total) will be observed by the nar- 
row field instruments during 7— ray emission phase. The 
UVOT is sensitive to magnitude 24 in a 1000 second ex- 
posure (For a linear increase of the sensitivity with the 
exposure time, that means a sensitivity of magnitude 19 
in a 10 second exposure). These two telescopes are suffi- 
cient to detect the very early optical emission predicted 
here. 

In this work, the problem has been treated under the 
ideal MHD limit. In fact, magnetic dissipation may play a 
role (e.g. Fan et al. 2004c). Thus our treatment is a sim- 
plification of the real situation, and further considerations 
are needed to fully depict the physics involved. 
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